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In this paper we estimate the influence of nonlinear gluon evolution in the production of prompt 
photons at the LHC pp collider. We assume the validity of collinear factorization and consider the 
EHKQS parton distributions, which are solutions of the GLR-MQ evolution equations and describe 
quite well the DESY ep HERA data, as input in our calculations. We flnd that both single and 
double photon production are enhanced for low-pr photons and central rapidities, while this effect is 
absent for the high-pr photons. The implications of this effect for the Quark-Gluon Plasma searches 
and for the QCD background to Higgs are also discussed. 
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Prompt photons at high energies has provided a direct 
probe of the dynamics of the strong interactions Jj] . Usu- 
ally photons are called prompt if they are coupled to the 
interacting quarks. In the framework of QCD the dom- 
inant production mechanism for the prompt photons at 
Tevatron and LHC colliders is the Compton scattering 
q + g ^ ^ + q (See, e.g. Refs. 0, HI). It is clear that the 
cross section of such process is sensitive to the gluon dis- 
tribution in the proton or nuclei 4]. In particular, LHC 
probes the gluon distribution at Xg « 10~^ |3|. Usually 
the parton distributions in a proton are described by the 
DGLAP evolution equations @ and the cross section is 
calculated using collinear factorization. More recently, 
the fcT-factorization, which takes into account effects of 
finite virtualities and transverse momenta of the incom- 
ing partons, has been applied to calculate the prompt 
photon production cross section [1, 0|- In this case, the 
cross section is given in terms of the unintegrated parton 
distributions, which are solution of the BFKL or CCFM 
evolution equations, for instance. However, at high en- 
ergies (small-x), these linear evolution equations predict 
a strong growth of the gluon distributions, which implies 
a large density of gluons in this regime and the violation 
of the unitarity. Consequently, new dynamical effects as- 
sociated to the unitarity corrections are expected to stop 
its further growth (For recent reviews see Ref. Q). 

About 24 years ago, Gribov, Levin and Ryskin (GLR) 
0, followed by Mueller and Qiu (MQ) [l^, performed a 
detailed study of the small- a; regime and argued that at 
large densities the physical processes of interaction and 
recombination of partons become important and should 
be considered in the QCD evolution, implying a nonlin- 
ear evolution equation. In the last two decades, the solu- 
tion and possible generalizations of the GLR-MQ equa- 
tion have been studied in great detail (See e.g. Ref. 
Q). Currently, one believe that the small- a; gluons in 



a hadron wavefunction should form a color glass conden- 
sate (CGC) which is described by an infinite hierarchy of 
coupled evolution equations for the correlators of Wilson 
lines [llj. In the absence of correlations, the first equa- 
tion in the Balitsky-JIMWLK hierarchy decouples and 
is then equivalent to the equation derived independently 
by Kovchegov [T2|. The resulting BK equation is almost 
equivalent to the GLR-MQ equation for the unintegrated 
gluon distribution, with the solutions presenting similar 
characteristics |13|. Experimentally, there are strong ev- 
idences of the nonlinear (saturation) effects at DESY- 
HERA. In particular, the DESY ep HERA data in the 
small- a; and low-Q^ region can be successfully described 
in terms of saturation models with the measured 
cross sections presentin g th e geometric scaling property 
[l5l|. Moreover, in Ref. [16l] the authors have shown that 
adding nonlinear corrections to the evolution, based on 
gluon recombination as derived in [9l.[l0j. can improve the 
overall leading order (LO) fits to the HERA DIS data. 
Its main results are that when nonlinear terms are in- 
cluded, the resulting gluon distribution is reduced with 
respect to the solution of the LO BFKL equation and the 
slowing of the evolution leads to an enhancement of 
the small- a; gluon distribution at < 10 GeV^ relative 
to the LO DGLAP gluon distributions. It implies simi- 
lar enhancements in the cross sections which are strongly 
dependent on the gluon distribution with respect those 
calculated using the solution of the LO DGLAP evolution 
equation. In the charm quark production has been 
estimated at LHC energies and a substantial enhance- 
ment was predicted. As the prompt photon production 
cross section is directly dependent on the behavior of the 
gluon distribution, one can expect a similar enhancement 
of the cross section. 

Our goal in this paper is estimate the magnitude of 
this enhancement and determine the kinematic region 
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FIG. 1: (color online) Rapidity dependence of the ratio Ry 
for two different choices of the hard scale Q and minimum 
transverse momentum p™'" . 
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FIG. 2: (color online) Transverse momentum dependence of 
the ration Rp^, for two different choices of the hard scale Q. 



where this effect should be more important. We con- 
sider prompt photon production via the QCD processes 
q + q ^ J + g and q + g 7 + g in order to inves- 
tigate by which amount nonhnear effects contribute to 
the hard process. Our main focus is the production of 
prompt photons in the kinematical region accessible at 
CERN LHC energy. Moreover, as a by product, we 
also estimate the contribution of these effects for the 
double photon production, which constitutes an impor- 
tant QCD background to the Higgs (h) production chain 

pp — !■ hX — !■ J"fX. 

Lets now present a brief review of the prompt photon 
production. Two types of processes contribute to the 
cross section: the so-called direct piece, where the photon 
is emitted via a point like coupling to a quark, and the 
fragmentation piece, in which the photon originates from 
the fragmentation of a final state parton. As the second 
component can be almost completely reduced by isolation 
criterion used in the experimental data analysis, we focus 
our study only in the direct component, which provide 
a clean probe of the hard scattering dynamics. In this 
case, we have that the prompt photon production cross 
section is given by [l| 



da. 
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where xt = 2pt/\/s, y and pT are the rapidity and trans- 
verse momentum of the produced photon, fi(x,Q^) are 
the parton densities, xi and X2 are the momentum frac- 
tions of the partons involved in the hard process. In this 
case we have that = ■^1^^'= " and a;',"™ — — 



Q is the hard scale and ^ are the partonic cross sec- 
tions, which are perturbatively calculable The con- 
tributing LO subprocesses are qg qj (Compton), 



qq gj (annihilation), followed by the subdominant di- 
agrams qq 77 (pure EM), gg 77 and gg — > gj. The 
correspondent LO matrix elements and partonic cross 
sections can be found in Refs. [11,113. Using Eq. © 
to calculate the cross section, we implicity are assuming 
the validity of the coUinear factorization. It is important 
to emphasize that factorization breaking isjjredicted by 
CGC physics in several processes (See e.g. [lj|). In what 
follows we calculate the prompt photon production cross 
section considering the EHKQS ^ and CTEQ6 ^ sets 
for the parton distributions and estimate the ratio among 
these predictions in order to quantify the magnitude of 
the enhancement expected at LHC. We analyze the ra- 
pidity and transverse momentum dependence of this ra- 
tio. As the EHKQS sets are only evolved to leading order 
(LO), we work the cross sections also at leading order. 
It implies that our predictions should be considered an 
upper bound, since the NLO small- a; gluon distributions 
are typically reduced relative to LO, implying a smaller 
enhancement. 

In Fig. [T] we present our estimates for the rapidity 
dependence of the ratio Ry defined by 
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We consider two different values for the factorization and 
renormalization scales: Q = pt and Q = 2pT, where pt 
is the photon transverse momentum. Moreover, in or- 
der to calculate the differential cross section da/dy we 
have assumed two distinct values for the minimum pho- 
ton transverse momentum: p™*" — 1 and 2 GeV. The 
results show an enhancement of low-px photons - the 
lower pjp™ the more significant enhancement of prompt 
photon production. We also notice a larger enhancement 
for central rapidities. The choice of the hard scale also af- 
fects the magnitude of the nonlinear effects, by choosing 
Q = Pt one has a bigger enhancement. In particular, for 
p-mm _ ^ QgV and Q = Pt, we predict a factor of two for 
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FIG. 3: (color online) Rapidity dependence of the ratio Ry 
calculated considering only the qq 77 and gg 77 sub- 
processes. 
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FIG. 4: (color online) Diphoton invariant mass dependence 
of the ratio Rai^ . 



the enhancement of the rapidity distribution associated 
to the nonhnear effects in the central rapidity region. 

In Fig. [2] we present our estimates for the transverse 
momentum dependence of the ratio Rpr^ defined by 
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Similarly to our previous analysis we also consider two 
different choices for the scale Q. Moreover, we have in- 
tegrated the rapidity distribution in the range \y\ < 3, 
which will be studied in LHC. Wc have that an enhance- 
ment is predicted in the region of small transverse mo- 
mentum pt < 10 GeV. At larger scales the nonlineari- 
ties die out since these terms are proportional to 
Consequently, the EHKQS gluons become similar to the 
CTEQ6 one, so that the enhancement disappears at large 
Pt- This behavior is also expected in charm production 
[17] . From the Fig. [5]we can see that the nonlinear effects 
imply an enhancement of 20 % in the small-pT region if 
we assume Q = 2pT. On the other hand, for Q = pt this 
enhancement can be larger than 100 %. The enhance- 
ment in the low range of transverse momenta has impor- 
tant implications, mainly if we remember that the magni- 
tude of prompt photon production in proton-proton col- 
lisions is used as baseline to estimate nuclear medium 
effects in nucleus- nucleus collisions |21j]. In particular, a 
signature of quark-gluon plasma formation is the thermal 
photon production, which manifests by an enhancement 
in the inclusive photon spectrum at px values below « 15 
GeV at LHC (For recent reviews see Refs. [22, 23]). The 
pre-requisite for the extraction of the thermal signal is a 
precise control of the photon production rate in proton- 
proton collisions, which can be strongly modified by the 
nonlinear effects, as demonstrated above. Consequently, 
our results indicate that a detailed study of the prompt 
photon production in pp collisions at LHC is necessary 
before to consider the thermal photon production as a 
precise signal of QGP formation. 



One can analyze the contribution of the nonlinear ef- 
fects for the diphoton production, characterized by the 
subprocesses qq 77 and gg — > 77, which have been 
included in our calculations. Although the gluon-gluon 
initiated process is not the dominant one, it contributes 
significantly for the 77 cross section. This may imply in 
bigger nonlinear effects, since both initial state small- 
X gluons should be amplified in the region < 10 
GeV^. In Fig. [3] we present our results for the rapid- 
ity dependence of the ratio calculated including only 
the qq and gg contributions. We consider the same two 
choices of factorization and renormalization scales as be- 
fore: Q — Pt and Q = '2pT, and we have assumed the 
value of p™*" — 1 GeV for the minimum photon trans- 
verse momentum. We can see that the enhancement due 
to nonlinear effects is more significant here than for single 
photon production - probably due the the relative impor- 
tance of the two gluon initiated contribution. In fact, if 
we had considered only the gg — > 77 and neglected the 
qq — > 77 contribution (just for theoretical study), the Ry 
ratio would increase to almost 5 (1.5) for the Q — pt 
[Q — 2pt) choice. Anyway, in the full LO contribution 
some of this increased enhancement remains. 

This result motivates a more detailed analysis of dipho- 
ton production, which is considered the main background 
for Higgs production via gluon fusion (See e.g. Ref. [2^1). 
In terms of the diphoton invariant mass M^^ and photon 
rapidities yi and the cross section can be written as 

Q 
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where xi = ^[e^'i -I- e?'^] and X2 = ^[e'^i + e-y^]. 
In particular, we have calculated the diphoton invari- 
ant mass distribution, da/dMjj, and estimated the ratio 
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Rm^^. defined by 
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In Fig. [3] we present our predictions for the behavior of 
this ratio in the kinematical range of M^^ < 50 GeV con- 
sidering two distinct values of the hard scale Q. We have 
that in this range an enhancement is predicted, mainly 
at small values of the factorization scale Q and diphoton 
invariant mass. On the other hand, for larger values of 
Mj~^ the nonlinear effects can be disregarded. Since the 
typical pt of photons needed to produce a Higgs would be 
around 50 GeV each, the nonlinear effects here considered 
will not be important for the QCD background to Higgs 
(For a discussion of nonlinear effects in Higgs production 
see e.g. Ref. {25*1). On the other hand, the enhance- 
ment in the low invariant mass region is predicted so the 
nonlinear contributions could be tested in this range. 
In summary, in this paper we have investigated the 



prompt photon production considering the coUinear fac- 
torization and the EHKQS parton distributions, which 
are solutions of the GLR-MQ evolution equations and de- 
scribe quite well the ep HERA data. Our results demon- 
strate that the nonlinear effects implies a large enhance- 
ment of the cross section for single and double photon 
production, which could be tested at LHC. In particular, 
the magnitude of the single photon production at small- 
Pt should be constrained before using thermal photon 
production as a signature of the QGP formation. On the 
other hand, the nonlinear effects can be disregarded in 
the kinematical range where double photon production 
is an important background for Higgs production. 
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